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Abstract 

Background: Pulmonary arterial hypertension (PAH) is a devastating disease with significant morbidity and 
mortality. At the macroscopic level, disease progression is observed as a complex interplay between mean 
pulmonary artery pressure, pulmonary vascular resistance, pulmonary vascular stiffness, arterial size, and flow. Wall 
shear stress (WSS) is known to mediate or be dependent on a number of these factors. Given that WSS is known to 
promote architectural vessel remodeling, it is imperative that the changes of this factor be quantified in the 
presence of PAH. 

Methods: In this study, we analyzed phase contrast imaging of the right pulmonary artery derived from 
cardiovascular magnetic resonance to quantify the local, temporal and circumferentially averaged WSS of a PAH 
population and a pediatric control population. In addition, information about flow and relative area change were 
derived. 

Results: Although the normotensive and PAH shear waveform exhibited a WSS profile which is uniform in 
magnitude and direction along the vessel circumference at systole, time-averaged WSS (2.2 ± 1.6 vs. 6.6 ±3.4 dynes/ 
cm 2 , P = 0.018) and systolic WSS (8.2 ± 5.0 v. 20.0 ± 9.1 dynes/cm 2 , P = 0.018) was significantly depressed in the PAH 
population as compared to the controls. BSA-indexed PA diameter was significantly larger in the PAH population 
(1 .5 ± 0.4 vs. 0.7 ± 0.1 cm/m 2 , P = 0.003). 

Conclusions: In the presence of preserved flow rates through a large PAH pulmonary artery, WSS is significantly 
decreased. This may have implications for proximal pulmonary artery remodeling and cellular function in the 
progression of PAH. 
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Background 

Pulmonary arterial hypertension (PAH) is a devastating 
disease with significant morbidity and mortality [1]. 
Classified as a subgroup of pulmonary hypertension by 
the World Health Organization, it includes persistent 
hypertension of the newborn, hypertension secondary to 
congenital heart disease, and idiopathic PAH [2]. For all 
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categories, survival rates are poor. Children with 
idiopathic PAH (comprising 50-60% of those with PAH) 
have demonstrated 5-year survival rate of 75% in both 
European studies as well as the US REVEAL (Registry to 
Evaluate Early and Long Term PAH Disease Manage- 
ment) study [3-5]. 

Research efforts are complicated by the multifactorial 
nature of the disease [6-14], involving maladaptive pul- 
monary wall remodeling and changes in the pulmonary 
arterial hemodynamic environment [2]. It is postulated 
that initial injury to endothelial cells in the distal 
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vasculature cause pulmonary artery medial hypertrophy, 
adventitial thickening, and neo-intimal lesions [15]. A 
resulting consequence is a progressive increase in 
pulmonary vascular resistance (PVR) and mean pulmon- 
ary artery pressure [15]. Moreover, an elevated mean 
pulmonary artery pressure is thought to distend the 
proximal arteries and increase pulmonary vascular 
stiffness [16,17]. These changes in artery compliance and 
size can ultimately affect the flow waveform and viscous 
hemodynamic forces at the artery walls. 

A key link in these events is the quantification of how 
wall shear stress (WSS) - the primary mechanical force 
affecting cell mechanotransduction - changes in disease 
conditions [18]. The WSS mechanism is attributed to 
the observation that high flow pulsatility promotes in- 
flammatory and proliferative cell expression [19,20]. 
While previous studies have quantified changes in mean 
pulmonary artery pressure, PVR, pulmonary vascular 
stiffness, arterial size, flow waveforms, and flow fields in 
the presence of PAH, WSS has not been evaluated quan- 
titatively in the pediatric PAH patient [6-14]. In addition, 
there is increasing evidence in the systemic vasculature 
that low WSS is a promoter of increased wall stiffness 
and atherogenic vascular states, and is an independent 
predictor of cardiovascular mortality [21,22]. Since WSS 
is reported to regulate transcriptional events in vascular 
remodeling, its quantification may further elucidate the 
complex etiology involved in PAH. In this study, we aim 
to quantify the local, temporal and circumferentially 
averaged WSS by cardiovascular magnetic resonance 
(CMR) in the right pulmonary artery (RPA) of a 
pediatric PAH population and a control population. In 
addition, we assess whether there is a difference in the 
measured WSS values between these groups. 

Methods 

With the approval of the Institutional Review Board at 
the Children's Hospital Colorado, all available CMR 
studies on patients with PAH at our institution were 
retrospectively analyzed. Patients with prior history of 
pulmonary artery surgery, right ventricular surgery, pul- 
monary arterial stenosis, pulmonary valve insufficiency, 
and chronic thrombolic pulmonary hypertension were 
excluded. In cases in which multiple studies were 
performed on a single patient, we chose to analyze the 
study with the best delineation of the RPA lumen and a 
high signal-to-noise ratio phase contrast data. Patients 
with a poor slice position through the RPA, sternal wire 
susceptibility artifacts, or poorly defined artery lumen 
boundaries were excluded (n = 5). For the normotensive 
group, we included patients who were referred for 
syncope with normal CMR, vascular rings, and mild left- 
sided cardiac anomalies that include non-stenotic bicus- 
pid aortic valve and uncomplicated coarctation repair. 



Demographic data for control versus PAH subjects are 
shown in Table 1. 

CMR Protocol 

A fast low-angle shot gradient echo sequence was used 
to obtain retrospectively gated tissue intensity and phase 
velocity maps encoded in the through-plane direction 
(1.5 T Siemens Magnetom Avanto). Double oblique im- 
aging slices were positioned between the proximal and 
the first branch of the RPA, orthogonal to the vessel 
long axis. A typical sequence used temporal resolution 
of 14-28 ms, echo times of 2.2-3.5 ms, and a flip angle 
of 25°. Depending on patient size and field of view, the 
cross sectional pixel resolution was 0.82 x 0.82 to 1.56 x 
1.56 mm with slice thicknesses of 5 mm. Velocity encod- 
ing values were adjusted according to the maximum vel- 
ocities encountered during scout sequences to optimize 
the velocity map resolution (typical values ranged from 
200 to 250 cm/s). 

Post-processing and RPA size 

The RPAs were carefully segmented (Figure 1) over the 
cardiac cycle using a semiautomatic level-set method 
(Segment, Medviso), and exported to a previously 
described Matlab (Mathworks, Inc., Natick, MA) pro- 
gram developed to record time-resolved RPA diameter 
measurements, blood flow (positive, negative and net), 
peak blood velocity (V max ), and the temporal and spatial 
WSS fields [23,24]. Pulmonary artery diameters were 



Table 1 Demographics and flow quantification of PAH 
study patients 





Control (n = 4) 


PAH (n = 25) 


P 


Age, (y) 


19± 12 


12 + 5 


0.3 


Sex, (M:F) 


2:2 


11:14 




BSA (m 2 ) 


1.7 ±0.7 


1 .2 ± 0.3 


0.14 


RPA diam (cm) 


1 .2 ± 0.5 


1 .8 ± 0.7 


0.19 


RPA/ BSA (cm/m 2 ) 


0.7 ±0.1 


1 .5 ± 0.4 


0.003 


RAC (%) 


32 + 18 


21 ±11 


0.15 


Vmean t avg (cm/s) 


30 + 12 


17± 11 


0.05 


Vmean systole (cm/s) 


83 ±23 


58 ±34 


0.07 


Qavg (L/min) 


2.6 ± 1 .2 


2.4 ± 1.0 


0.55 


Qsys (L/min) 


9.0 ±4.4 


8.9 ±3.7 


0.68 


V max (m/s) 


1 .3 ± 0.7 


0.8 ± 0.5 


0.14 


RF(%) 


0.9 ± 0.9 


1.1 ±2.0 


0.47 


WSS t avg (dyne/cm2) 


-6.6 ± 3.4 


-2.2 ± 1 .6 


0.018 



Values are mean ± SD; PAH = pulmonary arterial hypertension; BSA = body 
surface area (m 2 ); RPA = right pulmonary artery diameter (cm); RPA/BSA = 
normalized RPA diameter (cm/m 2 ); RAC = Relative Area Change (%); 
Vmean t avg = time-averaged mean velocity over vessel cross-section; 
Vmean syst0 | e = systolic mean velocity over vessel cross-section; Q avg = average 
flow rate (L/min); Q sys = peak systolic flow rate (L/min); RF = regurgitant 
fraction (%); V max = max velocity (m/s); WSS t _ avg = circumferentially averaged 
WSS over cardiac cycle. 



Truong et al. Journal of Cardiovascular Magnetic Resonance 2013, 15:81 
http://jcmr-online.eom/content/15/1/81 



Page 3 of 9 




Figure 1 Phase-contrast imaging through the right pulmonary artery. Segmented (a) intensity and (b) phase contrast images (RPA lumen 
shown in green). 



calculated by evaluating the ROI area at diastole and 
back calculating the effective diameter, assuming a circu- 
lar cross section (Figure 2). These were normalized by 
the body surface area, which was calculated by Haycock s 
method to obtain an indexed artery size for interpatient 
comparison [25]. The relative area change (RAC) of the 
arterial wall was calculated using the relation RAC = 
100*[(A max -A min )/A min ], where A max and A min represent 
the maximum and minimum measured cross-sectional 
areas. 

Flow & WSS quantification parameters 

The average flow rate (Q avg )> flow rate at systole (Q sys ), 
and regurgitant fraction (RF) were calculated according 
to convention. To facilitate the calculation of the shear 
fields and remove RPA bulk motion, the images were 



cropped to the RPA size and transformed from Cartesian 
to polar coordinates using oversampling (sampling at the 
vessel wall was set to double the original resolution) and 
bicubic interpolation to retain image integrity at the re- 
gion of interest (ROI) edges. As previously described, a 
linearized axial shear stress approximation was used to 
obtain the axial RPA shear field over all time-steps [23]. 
The axial WSS was then sampled at 45° increments 
along the RPA lumen to create a localized, position 
specific, waveform over the cardiac cycle. These posi- 
tions are referred to in the text and figures as S, A, I, P 
representing the superior, anterior, inferior, and posterior 
RPA lumen locations. WSS values averaged over the cir- 
cumference of the artery wall at flow systole ( WSS systo i e ) 
and averaged over the cardiac cycle (WSS t _ avg ) were also 
calculated as described previously [23]. 



Normotensive RPA 



Hypertensive RPA 



■ artery wall at systole 
wall deformation 
over cardiac cycle 




age = 14 
BSA= 1.33 m 2 
D = 2.00 cm 



age = 17 
BSA= 1.76 m 2 
D = 2.20 cm 



Figure 2 RPA area change and relative size. Normotensive (a) and PAH subjects (b) are shown with the calculated diameter (D), as back-calculated 
from diastolic cross-sectional area. 
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Measurement reliability 

In a subset of 10 random PAH and normotensive pa- 
tients, WSS t avg and WSS systo i e were calculated from 2 
separate segmentations of the RPA, to determine intra- 
observer variability (JD). Interobserver error was calcu- 
lated from RPA segmentations performed by two 
blinded researchers (AB and JD), and assessed with the 
values for WSS t _ avg and WSS systo i e . 

Statistical analysis 

All parameters are expressed as mean ± standard devi- 
ation. A Wilcoxon rank sum test was used to determine 
if parameter means were significantly different, where 
p < 0.05 was considered statistically significant. A regres- 
sion analysis similar to that previously reported [26] for 
healthy children was used to model the normal RPA 
diameter as a function of BSA (i.e. D = a*BSA 05 ) and 
95% confidence intervals were calculated to visualize 
deviation from normal values. 

Results 

Twenty- five PAH patients and 4 normotensive subjects 
with valid CMR images were analyzed. The mean age for 
the PAH group was 12 ± 5 years (ranging from 2 to 
20 years), with a median of 11.5 years. The mean age for 
the normotensive patient group was 19 ± 12 (ranging 
from 3 to 31 years). A summary of the demographics 
and bulk flow quantification parameters are shown 
Table 1. 

The normalized RPA size was significantly larger in 
the PAH population (p = 0.003) and evidence of a reduc- 
tion in RAC was observed (Figure 2), although not sta- 
tistically significant (p = 0.15). Regional WSS systo i e and 
the temporal WSS waveforms (Figures 3 and 4), were 
found to be markedly different. Figure 4 shows an 
example of the flow and regional WSS waveforms for a 
PAH patient and a normotensive subject. Both the 
normotensive and PAH shear waveform exhibit a WSS 
profile which is uniform in magnitude and direction 
along the vessel circumference at systole. However, the 
magnitude of the regional WSS systo i e differed drastically, 
as illustrated in the population means between the 
two cohorts for WSS t _ avg and WSS systo i e (Table 1 and 
Figure 5). 

Intraobserver variability analysis revealed that the per- 
centage error was 3.9% for WSS t _ avg and 5.1% for WSS sys . 
Interobserver variability was 2.7% for WSS t _ avg and 2.8% 
for VPSS^ 

Discussion 

There has been a paradigm shift away from the idea 
that pulmonary hypertension is merely a disease of the 
distal vessels [17,27]. Recent evidence has emerged 
warranting further examination of proximal pulmonary 




Figure 3 Regional wall shear stress in systole in the 
normotensive population (n = 4) compared to the PAH 
population (n = 25). P < 0.05 is indicated by '*'. 



artery remodeling in the disease pathway [17,27-29]. 
Proximal arterial dilation is a result of the second 
phase of proximal vascular stiffening in PAH, in which 
elastin function is lost, and the resulting dilated vessels 
must carry pressure load and maintain nearly constant 
diameter over the cardiac cycle (the geometric effects 
are observed via CMR in Figure 2) [17,30]. This in- 
creases RV load substantially and removes the dual- 
stage action of the right ventricle-pulmonary artery 
pump, requiring the RV to work even harder [31]. Such 
changes profoundly alter proximal flow dynamics, but 
it is uncertain how they affect overall pulmonary 
mechanotransduction as well as distal flow [17]. Alter- 
ation in pulsatile flow promotes inflammatory and pro- 
liferative cell expression [19,20]. Changes in WSS in 
the proximal arteries caused by changes in flow may be 
key to understanding the disease development. In this 
study, we are the first to describe depressed WSS 
values in a pediatric PAH population. 

Currently, the reference standard for diagnosing, 
monitoring, and evaluating therapeutic response in the 
presence of PAH is cardiac catheterization. The inva- 
sive nature of this procedure has led clinicians to look 
towards non-invasive imaging modalities for similar 
information. While echocardiography has multiple ad- 
vantages, CMR allows accurate assessment of ventricu- 
lar volume, mass, function, and evidence of myocardial 
fibrosis. CMR also has the potential to provide detailed 
WSS measurements, such as the case in atherosclerotic 
and aneurysm models [18,32,33], although this has not 



Truong et al. Journal of Cardiovascular Magnetic Resonance 2013, 15:81 
http://jcmr-online.eom/content/15/1/81 



Page 5 of 9 



Intensity Image 



Flow Rate vs Time 



Negative Flow 
Positive Flow 
Net Flow Rate 




Velocity Map 



Qsys =11.6 L/min 
Qavg = 3.33 L/min 



0.4 
Time [s] 

Local WSS (location shown by marker color in Vel. map) 



Peak Vel. = 125 cm/s 
Avg. Vel. = 96 cm/s 



-100 0 

cm/s 
Intensity Image 




Flow Rate vs Time 



— Negative Flow 
Positive Flow 
e— Net Flow Rate 




Velocity Map 



0.4 
Time [s] 

Local WSS (location shown by marker color in Vel. map) 




0.4 
Time [s] 

Figure 4 PAH patient morphology and WSS summary. Morphology, velocity, flow, and WSS summary for (a) an example normotensive 
subject and (b) an example PAH patient. The subjects are and BSA matched. Note the drastic overall WSS reduction in the PAH patient. 



been evaluated thoroughly in the pediatric patient. 
WSS, which is determined by the velocity gradient, 
cannot be measured by ultrasound Doppler modalities, 
and approximations using a Poiseuille (steady) profile 
have been shown to contain significant error. Previous 
efforts quantifying WSS in children have used compu- 
tational fluid dynamics to model the effect of PAH on 
WSS [34] and phase contrast (PC)-CMR has been used 
to quantify WSS in healthy and pre- and post-operative 
Fontan patients [35]. Tang et al. [36]. described WSS 



using computational fluid dynamics in PAH patients; 
however all but one was older than 18 years of age. 

We postulate that the observed proximal WSS de- 
crease is directly related to the increased size of the 
proximal pulmonary arteries in the presence of PAH. 
In preliminary studies, we used Color M-mode tissue 
Doppler echocardiography [27,31] to evaluate the differ- 
ences in RPA diameter in PAH children versus in 
normotensive children. Figure 6a depicts a combination 
of the previously collected RPA diameter data (n = 66 
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Figure 5 Cross-sectional average of l/l/SS sysfo/e in PAH patients 
as compared to normotensive controls ('*' indicates P < 0.05). 

P < 0.05 is indicated by '*'. 



PAH and n = 15 Control) and CMR data collected in this 
study (as plotted against BSA). The drastic departure of 
the PAH children from the 95% confidence bounds is 
evident with increasing BSA. An overlay of normal RPA 
diameters from Sluysman et al [26] (n = 496) reinforces 
this observation (Figure 6a, dashed line); additionally, 
Schiebler et al. [37] recently reported this enlargement 
trend in adult PAH. Given that compared to the normal 
subjects, PAH children had similar peak and average 
flow rates (Table 1) and that their RPA diameter was lar- 
ger, the result will be decreased WSS, as was observed 
for both WSS t _ aV g and WSS systo i e , Indeed, a strong trend 
toward decreasing WSS is seen as a function of the 
indexed RPA radius cubed (Figure 6b), which corre- 
sponds to the idealized pipe flow relationship between 



flow and diameter (i.e. WSS a [viscosity*flow/Diameter 3 ] 
[23]. These results suggest that dilatation is a progressive 
process in pediatric PAH and the consequence of 
maintaining equal flow through a larger artery is a de- 
crease in WSS. 

The results from our study are in concordance with 
previous efforts. Pelc et al. initially noted that main pul- 
monary artery flow at peak systole in PAH patients was 
more spatially heterogeneous and demonstrated a 
greater percentage of retrograde flow [38]. Building on 
this observation, a number of researchers have quanti- 
fied bulk flow hemodynamics in the pulmonary arteries. 
Among them, Ley performed a detailed study of 25 PAH 
patients and 25 volunteers to compare the peak velocity, 
average flow, time-to-peak velocity, velocity rise gradi- 
ent, and pulmonary distensibility [11]. In comparison to 
the volunteers, the PAH patients showed significantly 
reduced pulmonary velocities (P = 0.002), blood flow 
(P = 0.002), and pulmonary distensibility (P = 0.008). The 
patients also showed a shorter time-to-peak velocity 
(P < 0.001) with a steeper velocity rise gradient (P = 
0.002), which is consistent with higher pulsatility and 
increased higher harmonics of impedance. 

In addition, Morgan et al. has looked at the WSS 
forces in healthy pulmonary arteries [39]. As mentioned 
previously, WSS was found to be relatively constant 
throughout the cardiac cycle and axial, circumferential, 
and radial WSS averaged approximately 7 dynes/cm 2 . 
There was no statistical difference between the vessels 
for circumferential or radial WSS. Our study extends 
these previous investigations to determine the axial WSS 
forces in pediatric PAH patients. Table 1 summarizes the 
population demographics. Normalized RPA size, WSS tavg 
and WSS S y Sto i e were all found to differ significantly. The 




BSA, m 



— inverse cubic 
regression 
y=11.39*x 3 
r 2 =0.64, PO.001 



1 2 3 

Diameter/BSA (cm/m 2 ) 



Figure 6 Relationship between right pulmonary arterial diameter and WSS with BSA in controls and subjects with pulmonary arterial 
hypertension, (a) RPA diameter and WSS measurements in systole demonstrate a BSA dependence in both cohorts. The white line indicates the 
regression model (y = ax 05 ) for the control group, with the blue shaded region indicating the 95% confidence region; the dashed line indicates 
results (n = 496) from Sluysman et al. [26]. (b) WSS decreases rapidly as a function of the BSA-indexed RPA diameter. Solid markers indicate 
subjects with cross-sectional profiles plotted in Figure 2. The regression line in (b) is an inverse cubic regression, reflecting the inverse 
proportional relationship between WSS and the vessel radius cubed. Note: diameter measurements are augmented with data collected in 
previous echo-based studies [33]; WSS measurements are from CMR-only. 
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difference in the calculated WSS values between these two 
patient groups appears to be associated with artery dilation, 
as seen in the artery morphology comparison of matched 
patients in Figure 2 [11]. A case highlighting the effect 
artery size may have on WSS is illustrated by removing 
systolic flow as a variable between the BSA-matched 
normotensive and PAH patients shown in Figure 4. These 
patients have a similar systolic flow rate of 11.6 and 12.0 L/ 
min - with an indexed RPA size of 1.1 cm 2 /m 2 and 
2.7 cm 2 /m 2 , respectively. We postulate that as a result of 
this arterial size difference in the normotensive patient, and 
thus a change in spatial velocity gradients, the peak WSS 
magnitude is more than 4.5 times than that of the PAH 
patient. 

Gan used area distensibility and relative area change 
to predict mortality in PAH [7]. Besides finding that the 
cross-sectional area was much greater and the relative 
area change much less in PAH patients, relative area 
change was found to be the strongest predictor of mor- 
tality over distensibility, independent of pulse pressure. 
As well as measuring flow, pulse wave velocity, WSS, 
and distensibility in the proximal pulmonary arteries, a 
number of authors have observed highly inhomogeneous 
cross-sectional area flow profiles and retrograde flow in 
the main pulmonary artery of PAH patients [6,10,12]. 
Even though no correlation between retrograde flow and 
the degree of PAH has been found [10,12], these obser- 
vations suggest that PAH may influence some form of 
flow derangement in the right ventricular outflow tract 
[14]. Neither echocardiography nor 2D PC-CMR has 
shown any conclusive findings. However, Reiter et al. 
used time-resolved 3D velocity information to visualize this 
flow anomaly and found a vortex consistently formed in the 
main pulmonary artery of PAH patients, just below the 
RPA. It was found that the time persistence of this vortex 
was correlated to the degree of hypertension as measured 
by mean pulmonary artery pressure [14]. 

A number of concomitant factors contribute to changes 
in the pulmonary arterial pressure of PAH patients, in- 
cluding increased PVR, increased blood flow due to 
shunting, and alterations in pulmonary vascular stiffness. 
The influence of pathologic pulmonary pressures, PVR, 
and pulmonary vascular stiffness on proximal morpho- 
logic and hemodynamic factors may include pulmonary 
arterial dilatation, flow waveform alterations, and WSS al- 
terations [40]. Since WSS is reported to regulate transcrip- 
tional events in vascular remodeling, its quantification is 
important in order to model the complex etiology in- 
volved in PAH. PC-CMR is a promising tool for this appli- 
cation. It has been demonstrated to accurately and non- 
invasively quantify large artery spatial velocity fields, 
allowing for the measurement of spatial and temporal 
WSS [23]. Shear in PAH resembles an "atherosclerotic/ 
aneurysm phenotype" in that mean shear is greatly 



reduced proximally, potentially altering endothelial func- 
tion and mechanotransduction in the large vessels. Our 
group and others have shown that altered flow profiles 
(and thus, altered shear profiles) result in cellular ex- 
pression 

of inflammatory cytokines from more distal vessels 
[19,20]. This work is part of an ongoing effort to estab- 
lish a chain of events that eventually lead to worsening 
PAH; that is, dilation of vessels and distensibility lead 
to altered shear, which leads to altered cell expression, 
causing medial hypertrophy and endothelium prolifer- 
ation, elevating PVR, and ultimately culminating in 
right ventricular failure. 

Limitations 

Limitations to the study include those inherent to a 
retrospective study. We recognize that our CMR control 
population is small as this number is severely restricted 
by availability of RPA PC-CMR sequences in children 
without right heart disease. To address this, we included 
echocardiography results from healthy controls with no 
known heart disease, which shows significant differences 
in pulmonary arterial vasculature compared to pediatric 
PAH population. For the PC-CMR sequences, we in- 
cluded children with mild left-sided heart disease with 
no evidence of abnormal pulmonary artery pressure. 
However, we recognize the potential for mild pulmonary 
arterial flow abnormality caused by upstream anatomic 
anomalies. Efforts to enroll healthy children with no 
cardiac structural disease are ongoing. 

Conclusion 

We have demonstrated the use of PC-CMR to quantify 
hemodynamic WSS in PAH patients. Pulmonary artery 
dilation occurs at a significantly faster rate in PAH pa- 
tients, while systolic flow does not vary between popula- 
tions. Thus, the consequence of maintaining similar net 
flow through a larger artery is a WSS decrease in PAH 
patients. This may have implications for proximal pul- 
monary artery remodeling and cellular function. 
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